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ABSTRACT 



Recently Alard proposed to detect the shift of a star's image centroid, 8x, 
as a method to identify the lensed source among blended stars. Goldberg & 
Wozniak actually applied this method to the OGLE-1 database and found that 
7 out of 15 events showed significant centroid shifts of 5x ^ 0.2 arcsec. The 
amount of centroid shift has been estimated theoretically by Goldberg. However, 
he treated the problem in general and did not apply it to a particular survey 
or field, and thus based his estimates on simple toy model luminosity functions 
(i.e., power laws). In this paper, we construct the expected distribution of Sx for 
Galactic bulge events by using the precise stellar LF observed by Holtzman et 
al. using HST. Their LF is complete up to M/ ~ 9.0 (M v ~ 12), corresponding 
to faint M-type stars. In our analysis we find that regular blending cannot 
produce a large fraction of events with measurable centroid shifts. By contrast, 
a significant fraction of events would have measurable centroid shifts if they 
are affected by amplification-bias blending. Therefore, Goldberg & Wozniak's 
measurements of large centroid shifts for a large fraction of microlensing events 
confirms the prediction of Han and Alard that a large fraction of Galactic bulge 
events are affected by amplification-bias blending. 
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1. Introduction 

Experiments to detect Massive Astronomical Compact Halo Objects (MACHOs) by 
monitoring the light variations of stars undergoing gravitational microlensing events are 
being carried out by several groups (MACHO, Alcock et al. 1997a; EROS, Ansari et al. 1996; 
OGLE, Udalski et al. 1997; DUO, Alard & Guibert 1997). Since the lensing probability for 
a single source star is very low, these searches are being conducted toward very dense star 
fields such as the Large Magellanic Cloud and Galactic bulge. While searches towards these 
crowded fields result in an increased event rate, it also implies that many of the observed 
light curves include light from unresolved stars that are not being lensed: the blending 
problem. 

Depending on the source for the blended light, blending can be classified into several 
types. The first type, "regular blending", occurs when a bright source star registered on 
the template plate is lensed and its flux is blended with the light from numerous faint 
unresolved stars below the detection limit imposed by crowding. Regular blending affects 
the results of lensing experiments in various ways. First, it makes the measured event 
timescale shorter than the true one. Since the lens mass scales as M oc t%, the lens mass 
determined from the measured timescales will be underestimated and the lens population 
will be misinterpreted. In addition, since the optical depth is directly proportional to the 
summation of event timescales, i.e., r oc X^£e,i, the Galactic MACHO fraction determined 
from the optical depth without a proper blending correction is subject to great uncertainty 
(Di Stefano & Esin 1995). 

The second type of blending occurs if the source for the blended light is the lens itself: 
"lens blending" (Kamionkowski 1995; Buchalter, Kamionkowski, & Rich 1996; Buchalter 
& Kamionkowski 1997; Alard 1997). Besides the effects of regular blending, lens blending 
has several additional effects on the result of lensing experiments. First, because detecting 
events due to lenses close to the observer is comparatively more difficult than detecting 
events produced by lenses near the source, lens blending makes the optical depth depend 
on the geometry of the lens system. As a result, the matter distribution derived from the 
optical depth distribution deviates from its true one. Secondly, lens blending causes the 
lensing optical depth to depend on the lens mass function since more massive stars, which 
contribute more to the total optical depth, tend to be brighter, resulting in a larger blending 
effect (Nemiroff 1997; Han 1998). 

Finally, "amplification bias" blending occurs when one of several faint stars in the 
seeing disk below the detection limit is lensed, and the flux of the lensed star is associated 
with the flux from other stars in the integrated seeing disk (Nemiroff 1997). In current 
experiments, photometry is carried out by comparing template images, in which only very 
bright stars are resolved and registered, with a series of images taken of the same field. 
The result is that in amplification-biased events the brightest star appears to be the source 
because the lensed star is too faint to be resolved. To be detected, the amplification-biased 
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event must be highly amplified to overcome the high threshold flux from the brighter star. 
Therefore, the mean detection probability to detect events for each source star will be very 
low. However, if these faint stars comprise a significant fraction of the total number of 
stars, a considerable fraction of events might be amplification-biased (Han 1997b; Alard 
1997a). Moreover, due to the large amount of blended light from a bright star, the effects 
of blending for these events would be much more significant than those caused by other 
blending types. 

There have been several methods proposed to correct for the blending problem. The 
first method is to introduce an additional lensing parameter representing the residual flux 
from unresolved faint stars into the light curve fitting process. However, this method 
suffers from large uncertainties in the derived lensing parameters as a result of parameter 
degeneracies (Wozniak & Paczyhski 1997). Early-warning systems (for MACHO: Pratt 
et al. 1996; for OGLE: Udalski et al. 1994; for PLANET: Albrow et al. 1995) allow one 
to construct lensing light curves with high time resolution and small photometric errors, 
enabling one to detect small color shifts caused by blending (Buchalter et al. 1996). 
However, due to the narrow distribution of colors for Galactic bulge stars, the expected 
color shifts are small. Han (1997b) proposed to use the Hubble Space Telescope (HST) to 
provide blending corrections. By using the high resolving power of HST combined with 
color information from ground-based observations, one can identify the lensed source star 
in the blended seeing disk, thus the uncertainty in the derived timescale can be significantly 
reduced. However, this method requires costly HST times. One can also correct the 
blending effect statistically if the luminosity function(LF) of stars well below the detection 
limit can be constructed (Alcock et al. 1997b), but in this case we lose information about 
individual events. 

Recently Alard (1996) proposed to detect the shift of a star's image centroid, 5x, as 
a method to identify the lensed source among blended stars. Goldberg & Wozniak (1997) 
actually applied this method to the OGLE-1 database and found that 7 out of 15 events 
showed significant centroid shifts of 5x ^ 0.2 arcsec. The amount of centroid shift has 
been estimated theoretically by Goldberg (1997). However, he treated the problem in 
general and did not apply it to a particular survey or field, and thus based his estimates 
on simple toy model luminosity functions (i.e., power laws). In this paper, we construct 
the expected distribution of 5x for Galactic bulge events by using the precise stellar LF 
observed by Holtzman et al. (1998) using HST. Their LF is complete up to Mi ~ 9.0 
[My ~ 12), corresponding to faint M-type stars. In our analysis we find that regular 
blending cannot produce a large fraction of events with measurable centroid shifts. By 
contrast, a significant fraction of events would have measurable centroid shifts if they are 
affected by amplification-bias blending. Therefore, Goldberg & Wozniak's measurements of 
large centroid shifts for a large fraction of microlensing events confirms the prediction of 
Han (1997a) and Alard (1997) that a large fraction of Galactic bulge events are affected by 
amplification-bias blending. 
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2. Centroid Shift 

Consider a seeing disk in which many closely-spaced stars are located with positions 
and fluxes given by and F 0ti respectively. If one of these stars is gravitationally lensed 
and its flux is increased by A abs , the position of the center of light (CL) will shift toward 
the lensed star by an amount 

jU V /(A abs -1) + 1 A ohs ' [l) 

where A ^, s = 1 + /t(A a b s — 1) is the observed amplification, k = Fqj/ J2i Fa,i is the fractional 
light of the lensed source star before amplification, and (x) = (J2i^j -^o,i x i + -^ojXj)/ Hi Fo,i 
is the position of the CL before amplification (Goldberg 1997). Here the subscript "j" is 
used to designate the lensed star. 

In Figure 1, we illustrate the centroid shifts for a simple case in which two stars are 
blended and the light fraction of the lensed source star is k — 0.41. In the top panel, we 
present the point spread function (PSF) of each star. Since the separation between the two 
stars is small, the integrated PSF appears to be a single star (left middle panel) with its 
center at the CL (marked by a plus sign). In the rest of panels, we present the integrated 
PSFs and the position of CLs when A ahs = 3, 5, and 10 with respect to the original position 
of the CL when A abs = 1. As the amplification increases, one still cannot resolve individual 
PSFs, but the CL shifts increasingly towards the position of the lensed star. 



3. Expected Distribution of Centroid Shifts 

In this section we actually estimate the distribution of centroid shifts, f(Sx), which 
is expected for the microlensing experiment toward the Galactic bulge. Under current 
experimental strategy toward very crowded star fields, it is very likely that every event is 
affected by one of the three types of blending. However, the fraction of events affected by 
each type of blending is poorly known. Among the three types of blending, Han (1998) 
found that although a significant fraction of lenses might be stars, most of them are 
very faint, and thus lens blending has only a small effect on the result of microlensing 
experiments. We begin by estimating the individual distributions of centroid shifts caused 
by regular and amplification-bias blending, f reg (5x) and fa, mp (Sx), and then estimate f(Sx) 
for various relative combinations of both types of blending. 

To estimate f(5x) it is essential to know how densely stars are crowded and what their 
brightness distribution is, i.e., the precisely determined LF of Galactic bulge stars. We 
construct the LF by adopting the recent determination by Holtzman et al. (1998) using 
HST. Their LF is complete up to Mj ~ 9.0 (My ~ 12), corresponding to faint M-type 
stars. Although the expected contribution to the total event rate by stars fainter than this 
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limit would be very small, we extend the LF using the solar neighborhood LF from Gould, 
Bahcall, & Flynn (1997) as derived by HST imaging. We note that the faint part of both 
the Galactic bulge and disk LFs have similar shapes as seen by the good agreement between 
the two LFs in their overlapping region: 6.7 <^ Mj <^ 9.0. In Figure 2, we present the model 
LF. 

To simulate observations of microlensing events, we begin by using a Monte Carlo 
method to produce stars whose fluxes are distributed according to the model LF. At 
the same time, we assign positions to individual stars which are uniformly distributed 
throughout the seeing disk whose center lies at the CL. The average seeing in current 
microlensing experiments is ~ 2". However, when a star is located at the edge of the 
seeing disk, one can isolate its amplified image from that of the integrated seeing disk. We 
therefore set the size of the "effective seeing disk", the maximum unresolvable separation 
between stars, to 9 c g = 1".5. According to the model LF, there are on average ~ 10 stars 
within the effective seeing disk of angular area of 7r(# e ff/2) 2 . If the total flux of all stars 
within the effective seeing disk is greater than the detection limit imposed by crowding, we 
register this total as the template star brightness. Current observations reach detection 
limit when the total density of stars toward Baade's Window is ~ 10 6 stars/deg 2 . According 
to the model LF, this corresponds to M I ~ 3.1 mag. 

Next we select a lensed star. In the case of regular blending, the lensed star is the 
brightest star in the effective seeing disk. On the other hand, in the case of amplification- 
biased events, the lensed star is one of the relatively faint stars in the effective seeing disk. 
Once the source stars are chosen, we produce test events whose lensing impact parameters, 
f3, are distributed randomly in the range from to 1. Among these events, detectable events 
should satisfy the condition that the integrated flux at peak amplification is greater than 
the flux before amplification by a factor of 1.34, the threshold amplification, i.e., 



> 1-34, (2) 



0,i 



where A max = (f3 2 + 2)/(3((3 2 + 4) 1 / 2 . Finally, we compute the centroid shifts of individual 
detected events using equation (1). 

In the upper panel of Figure 3, we present the expected distribution of centroid shifts 
caused by regular (dotted line) and amplification-bias blending (dashed line). In the lower 
panel, we present the fraction of events with centroid shifts greater than a certain limiting 
value Sx\ im , i.e., 1 — J 5:Elim f(5x)dSx. In both panels, the distributions are normalized for a 
single source star, i.e., J °° f reg (Sx)dSx = J °° f amp (5x)d5x = 1.0 and the range < 5x < 1 is 
divided into 50 intervals. One finds several differences between the two distributions. First, 
while the peak of f Teg (Sx) occurs at 5x = and the distribution decreases rapidly with 
increasing values of 5x , the distribution f amp (Sx) peaks at 5x ~ 0'/2 and decreases gradually 
toward both smaller and larger values of 5x. In addition, while a significant fraction of 
amplification-biased events, ~ 50%, has centroid shifts greater than 5x\ im = 0.2 arcsec, the 
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fraction of regular blended events with 5x > 0.2 arcsec is just ~ 2%. 

Once the distribution of centroid shifts caused by each type of blending is obtained, 
the actual centroid shift distribution is obtained by 

f(8x) = (1 - J r amp)/reg(fe) + -^amp/amp (Sx) (3) 

with the known amplification-biased event fraction of ^ a m P - Han (1997a) estimated that 
J amp ~ 40% of events detected toward the Galactic bulge are affected by amplification-bias 
blending. However, to see the dependency of f(5x) on the value of .Tamp, we examine 
various amplification-biased event fractions, .Famp = 20%, 40%, and 80%. The resulting 
distributions of f(5x) and 1 — / falim f(Sx)dSx are presented in Figures 4 and 5, respectively. 
From these figures, one first finds that the distribution f(Sx) changes significantly 
for different values of J^amp- Secondly, in order to produce large centroid shifts for a 
significant number of events, a large fraction of Galactic bulge events must be affected by 
amplification-bias blending. Therefore, the recent finding by Goldberg & Wozniak (1997) 
that a considerable fraction of Galactic bulge events exhibit considerable centroid shifts 
confirms the predictions of Han (1997b) and Alard (1997a) that a substantial fraction of 
Galactic bulge events are affected by amplification-bias blending. 



4. Discussion 

The fact that events affected by amplification-bias blending produce large 5x while 
the centroid shifts caused by regular blending are small can be understood analytically 
in the following way. To produce centroid shifts large enough to be measured, events 
should satisfy two conditions. First, the event should be highly amplified. For a very low 
amplification event (A a b s ~ 1), the expected centroid shift will be small since — 1 ~ 
and thus f] ~ in equation (1). On the other hand, for a very high-amplification event 
(^4abs ~ oo), one finds 5x ~ |(x) — Xj| since the factor rj approaches unity. However, not all 
events that satisfy the first condition produce large centroid shifts. The second condition is 
that the lensed star should be one of the faint stars in the blended seeing disk. If the lensed 
star is the brightest one in the effective seeing disk and its flux dominates the flux over 
those from other faint blended stars, i.e., k ~ 1, the position of the CL before gravitational 
amplification will be very close to that of the lensed star, resulting in a small amount of 
shift (|(x) — Xj| ~ 0) since J2i^j -^o,i x i + ^oj^j ~ -^o,j x j an d J2i^o,i ~ F j. On the other 
hand, if the lensed source is very faint, i.e., n ~ 0, the light from the source star has 
negligible effect on the position of CL, resulting in high possibility of a large centroid shift. 
In summary, the expected centroid shifts for various extreme cases of amplification and 
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source star brightness are: 

rj ~ and 5x ~ for a low amplification event 

(^abs ~ 1) 

5x ~ | (x) — Xj | ~ for a high-amplification event 

with luminous source , , 

(A abs ~ 00, K ~ f) 

&c ~ | (x) — Xj I for a high-amplification event 

with faint source 

(A abs ~ 00, K ~ 0). 

For amplification-biased events, source stars are in general very faint (k ~ 0), mostly 
far below the detection limit. Despite their low luminosities, the fact that they are detected 
implies that the source stars are highly amplified (A^ ~ oo). Therefore, the conditions 
for large centroid shifts agree well with those for amplification-biased events. On the other 
hand, regular blended events do not meet these conditions. First, due to relatively small 
amount of blended light, regular blended events do not need to be highly amplified for 
detection (A ahs ~ 1). Although they can be highly amplified (A ahs ~ oo), the dominance of 
their fluxes (k ~ 1) over those from other faint blended stars will result in small Sx. 

As demonstrated by the large centroid shifts for a significant fraction of events, the 
effect of amplification-bias blending on the results of lensing experiments is important. 
However, the methods to correct for the blending effect mentioned in § 1 have various 
limitations in application. One very simple but very practical method to minimize 
the effect of amplification-bias blending is to monitor only very bright stars. With 
increasing reference image brightness, the required amplification for detection becomes 
higher, resulting in a lower probability of amplification-biased events. To show how the 
blending effect decreases with increasing reference image brightness, we simulate Galactic 
bulge events which are expected to be detected for various threshold reference image 
brightnesses. For each event, we compute the light fraction of the source star k and the 
timescale decrease factor r\. The distributions of k and rj are shown in the upper panel of 
Figure 6. For a given fraction of source star flux, the observed timescale is reduced by 

V = ^csAe — 2(1 — A~f n ) -1 / 2 — 2 ; A m - m = 0.34/k + 1. Here events are assumed to 
be detected as long as they can amplify the reference image flux by more than a factor 
of 1.34. However, highly blended events will have short t e s, resulting in low detectability. 
Therefore, we correct the distributions by the detection efficiency We assume that the 
detection efficiency is linearly proportional to the timescale decrease factor, i.e., e oc rj. 
The efficiency-corrected distributions /(/c) and j{r\) are presented in the middle panels. 
In the lower panels we present the distributions of the fraction of events with k > Ku m , 
1 — J Klim f(K)dK, and r\ > r)n m , 1 — Jq f(r)')dr)'. From these distributions one finds that 
under the current threshold reference image brightness of Mi ~ 3, the fraction of events 
with little blending effects (k ^ 0.9 or rj ^ 0.9) is ^ 10%. However, as the brightness of 
the threshold reference image increases, this fraction gradually increases until it becomes 
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~ 80% when only stars brighter than Mj ~ are monitored, which corresponds to the 
brightness of Galactic bulge clump giant stars (Paczyhski & Stanek 1998). The MACHO 
group (1997a) already applied this method and their optical depth determination is based 
on clump giant stars. However, by monitoring significantly fewer stars at a decreased event 
rate, the statistical precision of the lensing experiments will be lowered. A more general 
solution for the blending correction is provided by the rapidly progressing image subtraction 
technique (Alard & Lupton 1997b; Tomaney 1998) which is also being applied to detect 
microlensing events towards M31 by the Colombia- Vatt group (Crotts & Tomaney 1996; 
Tomaney & Crotts 1996). 

We would like to thank M. Everett and S. Gaudi for making useful comments. This 
work was supported by a grant 981-0203-010-1 from Korea Science and Engineering 
Foundation. 
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Figure 1: Illustration of the centroid shifts for a simple case in which two stars are blended 
and the fraction of light from the lensed source star is k = 0.41. In the top panel of the 
figure, we present the PSF of each star. Since the separation between the two stars is small, 
the integrated PSF appears to be a single star (left middle panel) with its center at the 
CL (marked by a plus sign). In the rest of panels, we present the integrated PSFs and the 
position of CL for A ahs = 3, 5, and 10 with respect to the original position of the CL before 
amplification. As the amplification increases, one still cannot resolve individual PSFs, but 
the CL shifts increasingly towards the position of the lensed star. 
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Figure 2: The model LF. We construct the LF up to Mj ~ 9 by adopting the recent 
determination by Holtzman et al. (1998) using HST, marked by filled squares, and further 
extend the LF by using the solar neighborhood LF from Gould, Bahcall, & Flynn (1997), 
marked by asterisks. We note that the faint part of the both Galactic bulge and disk LFs 
have similar shapes as seen by the good agreement between the two LFs in their overlapping 
region 6.7 £ Mj & 9.0. 
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Figure 3: Upper panel: The expected distributions of centroid shifts caused by regular 
[/reg(&0, dotted line] and amplification-bias blending [fa,m P (8x), dashed line]. Lower panel: 
the fraction of events with centroid shifts greater than a certain limiting value 8xn m , i.e., 
1 — J falim f(5x)d5x. In both panels, the distributions are normalized for a single source star, 
i.e., / °° f reg (5x)d5x = fa,m P (Sx)d8x = 1.0, and the range < 5x < 1 is divided into 50 
intervals. 
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Figure 4: The expected centroid shift distribution for events towards the Galactic 
bulge for various amplification-biased event fractions of JF amp = 20%, 40%, and 
80%. From the known distributions of centroid shifts due to regular, f[reg(8x)], and 
amplification-bias blending, /[amp(&r)], the actual centroid shift distribution is obtained 

by f(6x) = (1 - ^'amp)/reg(^) + ^amp/ampt^). 
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Figure 5: The distribution of the fraction of events with centroid shifts greater than a 
limiting value of 5x\\ ni for various amplification-biased event fractions of J^p = 0%, 20%, 
40%, and 80%. If all events are free from amplification-bias blending (jF amp = 0), only ~ 2% 
of events would have centroid shifts with 5x > 0.2 arcsec. On the other hand, to produce 
large centroid shifts for a significant number of events, a large fraction of Galactic bulge 
events should be affected by amplification-bias blending. 
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Figure 6: Correction for amplification-bias blending by monitoring bright source stars. 
Upper panels: The expected distributions of the fractional light k and timescale decrease 
factor 7] for various threshold brightness of reference image. Middle panels: same as upper 
panels, but the distributions are corrected by the detection efficiency. We assume that the 
detection efficiency is linearly proportional to the timescale decrease factor, i.e., e oc rj. 
Lower panels: The distributions of events with the fractional light n > Kn m (left panel) and 
the timescale decrease factor r\ > i]i im (right panel). 



